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Abstract. A sample of 97 galactic field Be stars were studied by taking into account the effects induced by the fast 
rotation on their fundamental parameters. All program stars were observed in the BCD spectrophotometric system 
in order to minimize the perturbations produced by the circumstellar environment on the spectral photospheric 
signatures. This is one of the first attempts at determining stellar masses and ages by simultaneously using model 
atmospheres and evolutionary tracks, both calculated for rotating objects. The stellar ages (r) normalized to 
the respective inferred time that each rotating star can spend in the main sequence phase (tms) reveal a mass- 
dependent trend. This trend shows that: a) there are Be stars spread over the whole interval ^ t/tms Ss 1 of 
the main sequence evolutionary phase; b) the distribution of points in the (t/tms, M/Mq) diagram indicates that 
in massive stars (M <; 12M©) the Be phenomenon is present at smaller t/tms age ratios than for less massive 
stars (M ^ 12M©). This distribution can be due to: i) higher mass-loss rates in massive objets, which can act to 
reduce the surface fast rotation; ii) circulation time scales to transport angular momentum from the core to the 
surface, which are longer the lower the stellar mass. 
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,£h 1- Introduction 



Correlations between the Balmer line emission width with 
Fsin i and the statistical tendency of Be- type emission line 
profiles to be present for low Vsmi yalues, while B e-shell 
type prevails at high Vsinj, inspired Istruvel l)l93l|) model 
of the Be phenomenon. This model underlies the most re- 
cent ones. The model also assumes that the re is a secularl y 
stable B-type stellar critical rigid rotator l|Tassoullll978|) . 
which builds an extended circumstellar envelope (CE) con- 
densed towards the equatorial plane by equatorial ejection 
of mass. Be stars are considered to be O, B and A spectral 
type non-supergiant stars that hav e shown at least onc e 
some emission in the Balmer lines ijjaschek et al.lli~98l|) . 
It has long been known that Be stars are fast rotators 
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and that they rotate at le ast 1.5 to 2 tim e s faster tha n 
B stars without emission jSl«t,t,fihalJ lT973 IZorecl r2004|L 
From a homogeneous set of Vsini para meters, although 
not co rrected for effects of fast rotation, IChauville et alJ 
(2001) concluded that Be stars rotate on average at an- 
gular velocity rates u) = Cl/Q c ~ 0.8. Stoccklcv (1968) 
pointed out that the Vsinz parameters can be systemat- 
ically underestimated if second order effects of fast rota- 
tion on absorption lines are neglected. Making allowance 
in the calculation of rotational line broadening for star dis- 
tortion and non-uniform surface temp erature in latitude 
l|von Zeipellll924albl) . Istoecklevl <|l968|) concluded that Be 
stars might actually be critical rotators. Thus, according 
to this author "...mild prominence activity or other minor 
disturbances lead to the ejection of matte r..." to fo r m the 
CE. These arguments were taken up by lOwocki I ll2004h 
and iTownsend et aD l)2004j) . In a study bv iFremat et alJ 
( 2005J) of rotational effects on fundamental stellar param- 
eters it is shown, however, that Be stars rotate on average 
atw~ 0.9. 
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ICrampin fc Hovlel l)l960|) suggested that the Be phe- 
nomenon occurs during the secondary contraction phase, 
where the surface rotation velocity has been spun up as 
a consequence of angular momentum conservation. Raw 
photometric color indices place Be star s near, or on, 
the TAMS (terminal-age main se quence) l|Schmidt-Kalerl 
ll964tlSchild fc Romanishinll976(l . However, this apparent 
location in the HR diagram of some Be stars is due in part 
to the continuum emission excess produced in CE and to 
the over-luminosity of the central objects carried by the 
rotationally-induced stellar geometrical defor mation and 
the concomitant gravitational darkening effect <lv on Zcipc 
1924albl:lsTettebak et alJl980tlMermillioJl982]:ISlettebak 



to another, which might suggest that the Be phenomenon 
can appear at any evolutionary stage during the MS ph ase. 



no- 



1985 t ICollins fc Sonnebornl Il977t ICoilins et all Il99ll 
Townsend et al.ll2004HFremat et alJr2005(l . 



If the Be phenomenon is an outgrowth of nearly criti- 
cal stellar rotation, one of the fundamental questions be- 
comes whether such a fast surface rotation is an innate 
property, or it is acquired at some stellar evolutionary 
phase. Two different phenomenological frames were put 
forward to tackle this question, which depend on whether 
Be stars are considered binaries or single stars. In bina- 
ries the Be phenomenon could arise after a Roche-lobe 
overflow event, when one of the components gains mass 
and angular momentum llPacketj F981: iHarmaned 1987; 
lGiesl200o]) . While this mechanism canno t entirely account 
for the observed frequency of Be stars l)Pols et alJll99li 
Ivan Bev er fc Va nbeverenlll997tK it can explain Be/X-ray 
binaries llCoel2000h . Therefore, let us assume that Be phe- 
nomenon concenrs only single stars only. In that case, the 
near critical equatorial velocity can either be an attribute 
of stars since their ZAMS (zero-age main sequence) phase, 
or a property that is acquired during their long-lived main 
sequence (MS) evolutionary phas e. Using moments of iner- 
tia of non-rotating stellar models . lllardorp fc Strittmatterl 
l|l970l) concluded that the initial rotation at ZAMS must 
be from 1 to 4% below the critical rate, for the star to 
become a critical rotator from core contraction in the MS. 
Since the fast rotation reduces co nsiderably the stellar mo- 
mentum of inertia, according to lEndall lll982h only 40% 
under-critical rotation at the ZAMS would be needed to 
accelerate the star to the critical rotation during the MS 
life span. For the star can reach critical rotation dur- 
ing the MS phase on the stellar surface also depends 
on the initial amount and internal distribution of an- 
gular momentum and on its loss and further redistribu- 
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Apart from the mentioned photometric results, few 
studies deal with th e ev olutionary status of Be stars . 
I.Iaschek et all l|l980h and iHubert-Delplace et all l|l982t) 
noticed a mild tendency of late Be type stars to belong to 
the giant luminosity class, while early Be type stars have a 
tendency to be dwarfs as if there w as a mass-age se l ection 
that underlies the Be phenomenon. IZorec fe Briotl (^997) 
found that the frequency of galactic Be stars against spec- 
tral type does not differ strongly from one luminosity class 



In a study of Be stars in open clusters . [Feinsteinl 
ticed an increase of the frequency of these objects near the 
middl e of the MS phase. Similarly, iFabregat fc Torre ionl 
( 200(1 concluded that the change of the frequency distri- 
bution of Be stars against the spectral type as a function 
of cluster age could be accounted for by assuming that 
the Be phenomenon occurs in the second half of the MS 
phase. 

Whil in the st ellar count of field Be stars by 
IZorec fe Briol l|l997j) the effects from CE-dependent over- 
luminosity and spectral changes due to fast rotation were 
taken into account, the binning of stars by luminosity class 
groups, which was meant to represent an evolutionary- 
dependent separation, cannot be justified en tirely for 
many stars of the sample. On the other hand, in lFeinsteinl 
l|l990j) and IFabregat fe Torreionl <|200C| . the sampling of 
Be stars in clusters against spectral type can be incom- 
plete. In these environments, Be stars are frequently de- 
tected photometrically. Since massive Be stars are scarce 
and the emission in the Balmer lines of stars cooler than 
B7 can be low, information on the appearance of the Be 
phenomenon in the relevant stellar masses may then be 
missing. Moreover, photometric and spectroscopic spec- 
tral types were not corrected for alterations due to the 
CE nor for changes induced by fast rotation. Thus, due 
to these shortcomings the conclusions drawn in those at- 
tempts are likely biased. 

Another way to undestand the evolutionary status of 
Be stars that may in principle solve the above inconve- 
niences, is to study in detail a statistically significant num- 
ber of individual Be stars, where in each star the perturba- 
tions produced by the circumstcllar emission/ absorptions 
on the observed spectra are considered and account is 
taken of the rotationally-induced effects. The aim of the 
present paper is thus to discuss the evolutionary stage of 
a well-observed sample of bright, field galactic Be stars, 
whose observational data were treated for all these de- 
viations. The fact that these objects are bright enough 
implies that their Be character is well recorded. This is of 
particular interest for those objects either among the more 
or the less massive stars of the sample, because in general 
they are not so numerous in cither of these extremes of 
mass. 



2. Observational data 

One of the main concerns related to the analysis of the 
observational data of Be stars is to remove circumstel- 
lar emission/absorption perturbations. We determine the 
fundamental paramet ers of these stars using B CD spec- 
trophotometric data ijChalonge fc Divanl ^352) . In this 
system the Balmer discontinuity is observed, which is 
characterized by two independent measurable quantities: 
the flux jump at A — 3700 A, D in dex and the mean spec- 
tral position of the discontinuity, Ai, presented in Ai— 3700 
A. The (Ai,Z?) parameters are strong functions of T e g 
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and log 3 l|Divan fc Zoredll982UZoredll986j) . It has been 
shown several times that for Be stars these two quanti- 
ties are free of circu mstellar extinction an d circumstellar 
emission/absorption ijZorec fc Briotlll99lj) . Be stars may 
show a Balmer discontinuity with two separated compo- 
nents. While the variable component due to the CE can 
be, in a given star, either in emission or in absorption de- 
pend i ng on the time they are observed ijMouitahid et alJ 
ll998Ul999l) . the constant component reflects the average 
photospheric prop erties of the observed stellar hemisphere 
i|Zorec et al1l2002() . 

The program stars are listed in Table ^ They were 
observed for more than 50 years at OHP (France) and 
ESQ (La Silla, Ch ile) with the Chalonge spectrograph 
llBaillet et al.lll973l) . a device specially conceived to ob- 
serve the stellar Balmer discontinuity. Since 1990 some of 
the program objects in the south hemisphere have been 
observed in CASLEO (San Juan, Argentina) with the 
Boiler & Chivens Cassegrain spectrograph, whose reso- 
lution in the low dispersio n mode is similar to th at of 
the Chalonge spectrog raph JCidale et alJl2000ll200lh and 
enables one to separate clearly the stellar from the cir- 
cumstellar Balmer discontinuity. The (Ai,D) parameters 
of the observed stars are listed in Table ^columns 2 and 
3). Their r.m.s deviations are on average ctd 0.005 dex 
and ^ 0.3 A. 

T he (Ai , D) were calibrated in the (T e g, log q) param- 
eters l|Divan fc Zoredll982llZoredll986j) . When these cali- 
brations are used for fast rotating stars, they must be con- 
sidered to represent the aspect angle-averaged properties 
of the stellar photosphere in the observed stellar hemi- 
sphere. Hereafter we call them apparent fundamental pa- 
rameters. The apparent (T e g, logg) sets are listed in Table 
n T c g and log g (columns 4 and 5). The Vsmi param- 
eters employed in the present work (column 6) are from 
IChauville et all l)200j) and Fremat et al. (2005). The listed 
Vsmi parameters wer e obtained using classical models 
of stellar atmospheres l|Stoecklev fc Mihaladll973|) . where 
the variation of the limb-darkening coefficient with fre- 
quency in the line is taken into account, so that they can 
be considered free of underestimations in duced by the use 
of co nstant limb-darkening coefficients ijCollins fc Truaxl 
19951) . Nevertheless, these Vsmi parameters must also be 
considered apparent. 

Fig. n shows the HD diagram of the observed Be stars 
given in terms of the observed BCD (Xi,D) parameters. 
In this diagram we can see the tendency mentioned by 
IJaschek et ail (|l980j) and lHubert-Delplace etaP l|l982|) of 
late type Be stars (cooler than B5) to be on average 
slightly more luminous than early Be stars. Although the 
(Xi,D) parameters of stars presented in Fig. ^ can be 
considered free of CE perturbations, they are "apparent" 
quantities because they still need to be treated for effects 
induced by fast rotation. Our stellar sample has a spectral 
type distribution that mirrors quite well th at of the whole 
know n Be star population near the Sun ijZorec fc Briotl 
119971) . The results we obtain with them can then be con- 
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Fig. 1. HR diagram of the program Be stars in terms of 
the BCD (Xi,D) parameters. Each curvilinear quadrilat- 
eral represents a MK spectral type-luminosity class group. 
The vertical strips demarcate the spectral types, which are 
given in the abscissa. The luminosity classes corresponding 
to the horizontal strips are marked in the quadrilaterals 
on the right side of the diagram 

sidered to represent fairly well the properties of this entire 
population. 

Table 1. Program stars, observed and derived parameters 
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3. Method 

3.1. Stellar atmosphere models of rotating stars 

The models of rapidly rotating stars used in the present 
work to describe the aspect angle average spectroscopic 
and spectrophoto metric charact e ristics of early type stars 
are described in iFremat et alJ l)2005|) (calculation code 
fastrot). They correspond to objects with overall rigid 
rotation and take into account their geometrical deforma- 
tion as described by equipotentials in the Roche approx- 
imation. Allowance is also made for changes of the polar 
radius and the bolometric luminosity produced in the stel- 
lar core. The lowering of the bolometric luminosity was 
considered to be related to the mass-compensation effect 
of rig idly rotating stellar cores l)Sackmannlll970t lOement] 
[1979J). The non- uniform effective temperature distribution 
with latitude follows the Ivon Zeipell il924albh theorem as 
far as high enough temperatures are concerned. For lo- 
cal effective temperatures lower than 8000 K we use d the 
gravitational darkening calculated bv lClaretl l|l998|) . The 
calculation code fastrot enables us to calculate spectral 
lines and energy distributions. We can then estimate the 
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Fig. 2. Model (Ai,£>) parameters against the rate fi/O c 
for different unperturbed pairs (T c ff, logg) and several as- 
pect angles i 

changes produced by rotation on the Balmer discontinuity 
and the (Ai, D) parameters that describe it. 

3.2. Relation between "apparent" and "parent 

non — rotating counterpart" stellar parameters 

The main purpose of the present paper is to infer stellar 
fundamental parameters that may give us some insight 
into the most plausible evolutionary state of the studied 
Be stars, once the observed quantities are treated for the 




Fig. 3. Functions Ft and Fq for t/tms = 0.5 and two 
stellar masses 



first order rotational effects. To derive the actual stellar 
mass M of a fast rotator, we assume that the observed set 
of parameters (Ai, D, Vsini) is affected by two types of 
rotational effects. There are direct changes related to the 
stellar geometrical deformation and the consequent non- 
uniform surface temperature and gravity distributions. 
Moreover, there are effects related to changes that the ro- 
tation produces on the evolution of stars. In order to take 
both types of effects into account, we proceed in two steps. 
Our models are built as a function of effective tempera- 
ture and surface gravity of homologous spherical stars: 
the same mass, but without rotation. So, in a first step, 
from the observed (Ai , D) quantities we derive the "parent 
non — rotating counterparts" (pnrc). In a second step, 
we use the pnrc (T e ff, logg) sets to derive the (T c ff,\og~g) 
quantities, which are the effective temperature and grav- 
ity averaged over the whole rotationally-deformed stellar 
surface. These average quantities are finally used as the 
entry parameters to the models of stellar evolution with 
rotation to infer masses and ages. For consistency with 
other works based on the use of fastrot, the nomencla- 
ture of fundamenta l parameters follows that adopted in 
iFremat et ail l|2005|) . 

The transition from apparent to pnrc parameters is 
carried out considering the following transformations: 



D = D (T eS ,\ogg) x F D (T cS ,\ogg, 
Ai - A°(T cff ,log 5 ) x F Al (T cffl log fl , 
Vsini = V" c (T cff ,log5) x 



fl e (T cff ,lo gg ^) 
R c (T eff ,log 9 ) 



(1) 



uj sin i 



where D and Ai are the observed BCD quantities, while 
D Q and A° are those of the rotationless homologous star; 
T e s and logg are the pnrc fundamental parameters; i is 
the stellar aspect angle; V c is the critical linear equato- 
rial velocity for rigid rotation; R c is the critical equatorial 
radius and R e the 'actual' equatorial radius at the rota- 
tional rate li — fi/f!, c . The factors Fr> and F\ 1 are func- 
tions calculated with fastrot for rigid rotators, which 
accounts for the geometrical deformation of stars, as well 
as for the corresponding aspect angle dependent external 
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non-uniform temperature and gravity distributions. Some 
examples of the behavior of D and Ai as a function of the 
angular velocity rate u> — Q/Q c and inclination angle i for 
several pnrc (T Q g,\ogg) pairs are shown in Fig El 

As in previous calculations of rotational effects on ob- 
servational quantities of B stars (|M^ dBrJg_Pextreniarm| 
197(ill972tlCollins fc Sonnebornll977tlCollins et alll99lt 



Townsend etaO 12004 IFremat et all 12005^ . the calcula- 
tions of (Ai,D) also show that rotational effects become 
conspicuous at rates u> 0.5. This limit is however lower 
for lower effective temperature. The (Ai,-D) parameters 
are double valued functions against uj at low enough pnrc 
effective temperatures (T e g & 13000 K). This occurs as a 
consequence of the gravitational darkening effect, which 
changes the average ionization balance in the stellar sur- 
face. As in stars with spectral types cooler than Al, the 
H absorption increases progressively over that of b-f 
transitions of neutral hydrogen, which carries a decrease 
of the value of D. Similarly, the Ai folding at high u> rates 
reflects the same behavior as in late B-type giant stars, in 
the sense that pressure broadening effects on Balmer lines 
makes Ai decrease more rapidly the lower the effective 
temperature as logg decreases. 

Since model evolutionary tracks are presented in terms 
of fundamental parameters averaged over the whole stel- 
lar surface (T e g , logg), a relation between these quantities 
and the pnrc (T e f{,logg) determined from has to be 
used to infer stellar masses and ages. Nevertheless, instead 
of using direct AQ = F(Xij , w) relations between pnrc and 
averaged X-quantities, we prefer to iterate relations like: 



Toff = T cff X F t (M,t,lj) 
log g = log g x F g (M,t,u) 



(2) 



which are essentially of geometrical nature and where the 
dependency on the stellar age r and mass M is small. Fig. 
shows the functions Ft and Fq for M = 3 and 20M© 
at t/tus = 0.5. 

The l^sin i parameter on the left-hand-side of the third 
relation in is intended to represent the true rotational 
parameter, i.e. the parameter corrected for underestima- 
tions induced by the gravity darkening l|Stoecklevl fl968: 
iTownsend et al.ll2004t IFremat et al.ll2005j) . Since this cor- 
rection depends on the pnrc (T e g,\ogg) and u>, it has to 
be iterated simultaneously as we search for the solution of 
the system 

To solve relations Q and J5J we need to specify either 
the angular v elocity rate u> or the inclination i. Since it 
was shown bv IFremat et al 1 l|2005h that most Be stars ro- 
tate at u> — 0.88([g4 06 , which is a distribution with a very 
low dispersion of angular velocities, we can adopt to = 
0.88 to solve the relations in the sought parameters M, r 
and i. In this paper we focus our discussion only on r and 
M. We also note that the solution of Q and J5J implies 
that we can translate (Ai, D) into (T c ff, \ogg). As empirical 
calibrations cannot be used, because of mixed rotational 
effects, we use model calculations. This implies that we 
do not have errors arrising from the procedure of trans- 
lating (Ai, D) into (T e ff, log g), but only with those passed 



from the observed BCD quantities onto the apparent fun- 
damental parameters. The propagation of empirical un- 
certainties in the determination of (t, M) is discussed in 
Sect.Q 



3.3. Evolutionary tracks of rotating stars 

Model tracks of stellar evolution with rotation calculated 
by Meynet & Maeder (2000, MM2000) for solar chemi- 
cal composition Z = 0.02 were done for different initial 
(or ZAMS) true equatorial rotation velocities V . In these 
models it is also assumed that in the ZAMS the stars start 
evolving as rigid rotators. The use of evolutionary tracks 
with rotation has two difficulties: 

a) no information exists on what V should be adopted 
to interpolate stellar masses and ages. IZorec et alJ 1 20041) 
have shown that the true equatorial velocities of dwarf Be 
stars have a quite flat distribution against spectral type 



around V ~ 300 km s" 1 (V ~ 350 km s^at BV0 and 
V ~ 270 km s _1 for BV9). Calculations of internal an- 
gular momentum redistribution foresees that in the first 
1 to 2% of the MS lifetime an initial flat internal angu- 
lar velocity distribution transforms into a step-like one, 
where, depending on the mass, the rotation in the stel- 
lar c ore becomes 20% to 40% faster than in the enve - 
lope llPenissenkov et al J IToQ^ iMevnet fc Maederl EoOfi. 
Since classical Be stars have masses that range from 3 
to 30Mq, we should then use models of stellar evolution 
calculated for somewhat higher mass-dependent initial ve- 
locities, ranging from V — 340 at 3M to 420 km s _1 for 
30M Q ; 

b) there is some evidence for internal angular momen- 
tum redistribution in the pre-main-sequence (PMS) evo- 
lutionary stages o f stars with masses from 0.1 to IOMq 
I Wolff et alJ 120041) which implies that stars can start 
evolving from the ZAMS as differential rotators. If so, and 
depending on the PMS evolutionary characteristics of in- 
dividual stars, the amount of rotational energy stored by 
them could be hig her than the li mit imposed by the crit- 
ical rigid rotation l|Tassoullll978|) . In such a case there is 
a much higher mass-compensation effect on the core bolo- 
metric luminosity and there may be more consequences on 
the stellar evolution than those accounted for in models 
used in the present work. 

Facing the quoted unknowns on the internal rotation 
of the studied stars at the ZAMS and on their initial equa- 
torial velocity, we estimate the effect on the mass and age 
estimates in the physical framework defined by the exist- 
ing calculations of stellar evolution with rotation. 

First, we obtained stellar ages and masses with the 
evolutionary tracks without rotation for Z = 0.02 
l|Schaller et al.lfl992f) . Then, we derived the same quan- 
tities using the evolutionary tracks for V — 300 km s _1 
and Z = 0.02 of MM2000. These last correspond to initial 
angular velocity rates that range from ui = 0.79 at Al = 
3M Q to uj = 0.52 for M = 30M Q . They produce smaller 
effects on the mass and age estimates than the slightly 
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higher velocities V = f(M) needed to produce the aver- 
age main sequence V ~ 300 km s" 1 after having undergone 
rapid initial internal angular redistribution in the ZAMS. 
The enhanced values V a = f(M) imply rotational effects 
on stellar evolution scaled in terms of initial rates ranging 
from lu = 0.87 at M = 3M© to 0.70 for 3OM . This dif- 
ference in the initial rates w may have some significance, 
since rotationally-induced effects on the stellar fundamen- 
tal parameters increase rapidly once w >0.8 and th ey are 
stronger for lower stellar mass l)Fremat et alJEof)^ . 

In order to estimate an order of magnitude of the ef- 
fects on the estimates of mass and age by initial equatorial 
velocities larger than those used in the published models, 
we re-scaled the existing evolutionary tracks as a function 
of the abovementioned mass-dependent values V = f{M). 

3.4. Re-scaled evolutionary tracks 

The calculations carried out bv lEndal fc Sofial l)l979j) and 

more recently by 

iHeger fc Larigerl ll2000l) : IMevnet fc Maederl i200Cl l2002h 
and iMae^eTfcMevnetl l|200lj) show that rotation intro- 
duces several changes in the evolutionary tracks com- 
pared to those for non-rotating stars. The characteristics 
of these changes, under the assumption that the stars start 
evolving from the ZAMS as rigid rotators, depend on: 
the adopted mass-loss rates, the initial conditions such as 
chemical composition and the initial rotational velocity V 
and on the mechanisms of angular momentum redistribu- 
tion inside the star. In Be stars, the average mass-loss rate 
Ss 10~ 9 M Q yr , which encompasses winds and discrete 
mass ejections, cannot lead to sensitive deviations from 
the evolution with the time-dependent variation of stel- 
lar mass already foreseen in the existing calculations. In 
this work, those changes of mass are assumed to be the 
same as that calculated for objects evolving with V = 300 
km s . However, depending on the initial value of the ro- 
tational velocity and the further phenomenon of angular 
momentum redistribution, there are at least three other 
outstanding changes in the evolutionary tracks of rotat- 
ing stars in the MS phase that interest our fundamental 
parameter determination: 

i) In the (logL/L , T ff)-plane the tracks are slightly 
shifted and rotated, so that for a given mass the starting 
point in the ZAMS is located at a lower temperature and 
luminosity, which re veals the rotatio nally-induced mass- 
compensation effect l|SackmanrJI 1 9 70|) : 

ii) The MS phase is prolonged to higher luminosities 
than in the non-rotation models due to the enlargement 
of the H-content in the convective core, which is produced 
by t he mixing processes that fuel it with fresh hy dro- 
gen l|Heger fc Langerll2000t IMevnet fc Maederl |2000|) . On 
the other hand, as evolution proceeds in the MS phase, 
in rotating stars there is a more sensitive change of the 
mom e nt of inertia than in no n-rotating objects (Endal 
Il982t IMevnet fc Maederl |2000|) . which leads to an en- 
hanced stretching of the star. The MS phase can then end 



5.5 




log T eff 

Fig. 4. Evolutionary tracks for different initial velocities 
V . 'dotted 1 lines are for evolutionary tracks with V D = 
km s _1 ; l fuW lines are for V = 300 km s _1 ; 'dashed' 
lines correspond to V = f(M) so that V — 340 km s _1 
at 3M n and 420 km s" 1 for 30M ra 



up at lower surface-averaged effective temperatures than 
in non-rotating models; 

Hi) The overall evolutionary MS life span outlasts the 
non-rotating case, on the one hand because the levita- 
tion effect produced by the rotation m akes stars behave 
as is they had a l o wer effective mass llSackmanrJ ll970: 
lRodenheimerlll97lt IClemend Il979l iMoss fc Smith! Ii 982ft 
and on the other hand, because the mixin g processes cited 
in ii) increase the core-burnin g lifetime l|Heger fc Langerl 
l2000( IMevnet fc MaederlhoOOl) . 

These effects give a first order insight into the changes 
induced by the rotation, which actually encompass a se- 
ries of subtle global and local mechanisms of mixing and 
angular momentum redistribution, whose incidence on the 
stellar structure can be appreciated only through detailed 
calculations. In this work we use the final settings of such 
calculations carried out by MM2000 to re-scale the evo- 
lutionary tracks for initial velocities V a at will. We use 
the term 're-scale' instead of 'interpolate', because some 
effects calculated in detail have been published explicitly 
only for a restricted number of stellar masses and some- 
times they were typified by a single case that we had to 
generalize. Only the MS phase was treated in this way. 
The (logi/L Q , logTcff) evolutionary paths in the post- 
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MS part were constructed by translating the predictions 
given by model-evolution without rotation. Since by defi- 
nition Be stars are non-supergiant, the supergiant region 
will rarely be used in the present work. The results of this 
re-scaling operation is shown in Fig. 0] We see there: 1) 
evolutionary tracks for non-rotating stars (dotted lines); 2) 
model tracks by MM2000 calculated with V Q = 300 km s _1 
in all masses (full lines); 3) re-scaled tracks for a linear en- 
hancement of the initial velocity V above 300 km s" 1 as 
a function of mass given by V a = 300+ [3 x (M/M Q ) + 30] 
km s _1 (dashed lines). This relation is intended to give 
the initial velocities so that after the first « 10 4 yr in the 
MS phase the Be star true equatorial velocities become 
V ~ 300 km s- 1 as observed (cf. SectEHf. 

Changes of initial velocities from V = 300 km s _1 to 
V = 300+ [3 x (M/Mq) + 30] km s" 1 seem not to produce 
huge effects on the evolution path shapes. This saturation 
effect noted for high V values was already commented on 
by MM2000 for hot rotating stars, where increased mass- 
loss rates due to a faster rotation partially suppress the 
effects related to an enhanced outward transport of an- 
gular momentum. Our 'transcription' of the phenomenon 
in all masses for whatever V a may perhaps be insufficient 
as it depends on a simple interpolati on among the few 
cases given explicitly in the literature l(Talon et alJll997t 
iMevnet fc Maedeill2000h . However, we note in Fig. H that: 
1) the breadth of the MS phase of rotating models is 
slightly enlarged for all OB stellar masses as compared to 
the rotationless case; 2) the stellar mass that corresponds 
to a given (logL/ Lq, logT c ff) point inferred in the lower 
half of the MS phase of rotating objects will be slightly 
higher than the same obtained with tracks of non-rotating 
models; 3) stellar masses inferred in the upper half of the 
MS will not depend sensitively on the type of evolution 
model used, if masses are lower than ~ WMq, while for 
stars with M ^ 10M Q the masses obtained from rotation- 
dependent tracks will be lower compared to those inter- 
polated with rotationless evolutionary paths. 



3.5. Stellar ages 

The evolution time scales depend on the total angular 
momentum and its progressive internal redistribution. Let 
us call r^g the time spent in the MS phase by a non- 
rotating star (from ZAMS to TAMS) and use the notation 
tms for the MS life span of its homologous rotating object. 
Assuming that stars have rigid rotation in the ZAMS, so 
that their initial equatorial velocity is V , MM2000 have 
shown that between tms & n d Tm S the following relation 
holds: 



tms 



= 1.0 + aV MS (M,V o ) 



(3) 



'MS 



where a = 0.0013 and Vms(M, V a ) in km s _1 is the surface 
rotational velocity averaged over the whole tms period. 
On account of the simplicity of relation |HjJ) and that ac- 
tually we do not have another way to scale stellar ages on 
evolutionary tracks for whatever initial values of V , under 



Table 2. Ratios 7 as a function of mass and the initial 
angular velocity rate iv = f(V /V c ) 



M/Mq 


U!o = 0.5 


0.8 


0.9 


20.0 


-0.022 


-0.008 


-0.008 


15.0 


-0.011 


0.021 


0.029 


10.0 


0.003 


0.056 


0.074 


5.0 


0.022 


0.113 


0.154 


3.0 


0.036 


0.170 


0.242 



the assumption that stars are rigid rotators in the ZAMS, 
let us write a similar relation to Q for the age t(oj) of 
a rotating object at any moment of its MS evolutionary 
phase and the respective 'rest' age r° = r(w=0): 



= 1.0 



V(t)dt 



(4) 



where V(t) is the time dependent surface velocity. After 
a rapid decrease that lasts 1 to 2% of the MS lifetime, in 
massive stars with mass-loss, V(t) shows a roughly 'linear' 
decrease with time (MM2000), while in stars with masses 
M < 10M Q , wher e mass-loss is negligible, it follows a 
'parabolic' increase l)Endalll982HMevnet fc Maederl2 0001 . 
We can then derive the following relation between the age 
fraction spent in the MS phase by a rotating star and the 
homologous non-rotating object: 



iH 

TMS 1 + l(r°/T MS ) 

in which wc have : 
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Vms(M,V ) 

Vo 



q(M,V ) 



aV a 



1 + q(M, Vo)aV 



(5) 



(6) 



where q(M,V ) = Vi/V < 1; Vi represents the surface 
equatorial velocity after the very initial short-lasting in- 
ternal angular momentum redistribution. Using the values 
of Vi/V and Vms(M, V n ) as a function of M and V from 
MM2000 and iDenissenkov et alJ (|l999^) we obtained the 
estimates of 7 that are given in Table El We see there that 
in many cases it is I7I < 0.1 so that to a good approxima- 
tion: 



(7) 



r(w) — tms ( 7 Jo 



MS 



On the other hand, the small values of |7| warrant that 
a very detailed representation of the function V(t) is not 
relevant to evaluate the age ratios JSJ). 

Since we will be dealing with uj ^ 0.8, most Be stars 
studied will have t(cj)/tms r °/ r MS- general devia- 
tions from (JJJ are small and, as expected, they are higher 
the lower the mass. 

For interpolation of stellar ages, we divided the 
ZAMS-TAMS time interval of each evolutionary path 
into 100 parts: 73, % =0,1,. ..100, so that tq =tzams an d 



8 



J. Zorec et al.: Evolutionary status of Be stars 



Tioo =ttams =TmS- Times Tj for the same i in tracks 
of two consecutive masses were considered representing 
stars in homologous evolutionary stages. Interpolations of 
masses and ages were then done in the (log L, logT ef f ) dia- 
grams according to this criterion of homologous evolution. 

A few studied objects fall in the HR-diagram strip of 
the secondary contraction phase. This does not allow us 
to decide if these objects are still in the MS phase, in the 
secondary contraction region, or of they are already in the 
giant branch. We treated them as if they still were in the 
MS, since the short time scales involved in the remaining 
phases do not change much the estimate of the age ratios 
t/tms- A small number of stars remain in the post-MS 
phase, even after all corrections of parameters for rota- 
tional effects. In the giant phase the evolution times were 
re-scaled using a relation similar to Q. 

4. Results 

4.1. Preliminaries 

The major motivation of the present paper is to infer the 
'present-day' age of stars already displaying the Be phe- 
nomenon. We cannot say if the phenomenon has already 
been present in a given star for some time, or if the star 
will display it up to the end of its MS phase. 

As alredy noted, Fremat et al. (2005) have shown that 
most Be stars rotate at u> ~ 0.88. Since this rate applies 
to stars that can be at different evolutionary stages in the 
MS phase, only those tracks that imply a state of the sur- 
face velocity that fit the condition u> ~ 0.9 at the required 
location of the star in the HR-diagram would be suitable 
to infer its mass and age. However, we do not know the 
individual initial velocities V to build the required model 
tracks. This can be solved partially by iterating the stel- 
lar mass and its V . Nonetheless, the operation requires 
a number of subtleties that are beyond the scope of the 
present work and will be developed elsewhere. The results 
obtained in this section will show a posteriori that adopt- 
ing an appropriate V for all stars is an approximation 
that suffices for the purposes of the present work. So, in 
the present paper we calculate masses and ages adopting 
the models by MM2000 for V = 300 km s -1 . We esti- 
mate the magnitude of possible uncertainties caused by 
the lack of knowledge of the specific value of V . To this 
end, we determine the (t/tmSi M/Mq) parameters of sev- 
eral HR 'test' points in the upper and lower half of the 
MS phase (chosen at hoc), as a function of evolutionary 
tracks dependent on different values of V . We assume 
that the 'test' points correspond to stars rotating at u) = 
0.88. A given set of test apparent parameters produce, as 
expected, inclination angle-dependent series of pnrc and 
averaged fundamental parameters. However, the relative 
changes from the use of various evolutionary tracks are the 
same. So, we assume that the test objects are seen at i ~ 
52°, or an average inclination of rotation axes oriented at 
random (52° ~ arcsin[sini = 7r/4]). For this specific angle, 
apparent X(i) fundamental parameters and the respective 



Table 3. Comparison of masses, ages and age ratios de- 
rived from evolutionary tracks without and with rotation 



N° 


T c ff 


log g 


M/Mq 


T (age) 


t/tms 




K 


dex 




years 




1 


Apparent parameters 


Evolution without rotation 




36324 


4.043 


24.78 


2.24x10 s 


0.350 




29119 


3.449 


24.03 


5.70xl0 6 


0.900 




27102 


4.100 


11.98 


5.53xl0 6 


0.350 




23529 


3.630 


11.93 


1.41xl0 7 


0.900 




14302 


4.157 


4.00 


5.72xl0 7 


0.350 




12169 


3.695 


4.00 


1.46 x 10 s 


0.900 


o 
z 


pnrc parameters for 


Evolution without rotation 




n/n c = o 


88 










38903 


4.185 


27.67 


2.68 xlO 5 


0.046 




30941 


3.563 


25.10 


5.16x10 s 


0.841 




29169 


4.234 


13.15 


4.23 xlO 5 


0.031 




25291 


3.751 


12.67 


1.16xl0 7 


0.807 




15484 


4.280 


4.28 


7.37x10 s 


0.054 




13176 


3.808 


4.19 


1.18xl0 8 


0.818 


q 
o 


Averaged parameters 


Evolution without rotation 




for Q/Q c ~- 


= 0.88 










36529 


4.037 


25.35 


2.24x10 s 


0.359 




29026 


3.437 


24.24 


5.67x10 s 


0.900 




27077 


4.089 


12.12 


6.29 x10 s 


0.407 




23478 


3.622 


12.15 


1.39xl0 7 


0.910 




14294 


4.149 


4.06 


6.62xl0 7 


0.425 




12163 


3.690 


4.09 


1.41 x 10 s 


0.922 


A 

4 


Averaged parameters 


Evolution with rotation 




for 0,/Q.c - 


= 0.88 


V = 300 km/s V masses 




36545 


4.046 


25.88 


2.31x10 s 


0.299 




28991 


3.411 


22.83 


7.60 x10 s 


0.893 




27089 


4.106 


12.61 


5.08 x10 s 


0.284 




23473 


3.614 


11.92 


1.64xl0 7 


0.863 




14294 


4.171 


4.27 


4.10xl0 7 


0.241 




12163 


3.690 


4.10 


1.59xl0 8 


0.842 


c: 


Averaged parameters 


Evolution with rotation 




for 0,/Q.c - 


= 0.88 


V = f(M) > 300 km/s 




36566 


4.058 


26.61 


2.00 x10 s 


0.253 




28989 


3.410 


22.75 


8.02x10 s 


0.898 




27094 


4.115 


12.85 


4.38 x10 s 


0.243 




23473 


3.614 


11.92 


1.69 xlO 7 


0.867 




14294 


4.176 


4.32 


3.44xl0 7 


0.204 




12163 


3.690 


4.10 


1.63 xlO 8 


0.846 



N° = block identifier to indicate in the text the type of 
evolutionary model used 



surface averaged X obey: X(i ~ 52°) ~ X (the Xs stand 
for logi, log g, logT c ff, etc.). 

The results thus obtained are displayed in Tabled The 
(T e ff,logg) parameters given in columns 1 and 2 of the 
1st block in Table [3] represent the observed, i.e. apparent, 
fundamental parameters which need to be treated for ro- 
tational effects. In columns 3 to 5 of the 1st block are 
given the masses, ages and fractions of age spent in the 
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MS as reflected by models of stellar evolution without ro- 
tation (Schal ler et al.ll992|) . In columns 1 and 2 of the 2nd 
block are displayed the respective sets of pnrc (T c g , logg) 
parameters of test points. They correspond to parameters 
the stars would have at rest. In columns 3 to 4 of the 2nd 
block we display the 'fictitious' quantities if the pnrc fun- 
damental parameters were used to derive masses and ages 
from evolutionary models without rotation. The pnrc pa- 
rameters (T ff, log g) given in columns 1 and 2 of t he 2nd 
block were derived using the iFremat et al] l|2f)f)5l) model 
atmospheres for rotating stars. Since evolutionary tracks 
of rotating stars are given in terms of fundamental pa- 
rameters averaged over the rotationally-deformed stellar 
surface, in columns 1 and 2 of the 3rd, Ath and 5th blocks 
we give the surface averaged effective temperatures and 
gravities of the test stars rotating at f2/f2 c ~ 0.88. In 
columns 3 to 5 of the 3rd block are given the masses, ages 
and fractions of MS ages derived using evolutionary tracks 
without rotation, while in the same columns of the Ath and 
5th blocks we give the parameters inferred using the orig- 
inal models by MM2000 with V Q = 300 km s" 1 in the 
ZAMS and the re-scaled evolutionary tracks for the V — 
f(M) meant to account for the average V = 300 km s _1 
of dwarf Be stars after the initial fast redistribution of the 
internal angular momentum (Sect. EH ■ 



The figures in Table |3] reveal that for mass estimates, 
the uncertainties from possible mismatches between V a 
and fi/fi c at the required location of the star in the HR 
diagram are not higher than 1M© for M il 10M© and 
they are much smaller for masses M <J IOMq. There may 
be, however, strong differences in the absolute age esti- 
mates. These differences have to be taken into account 
when comparisons must be done with ages of stars in en- 
vironments like clusters that were inferred from non- or 
slowly- rotating stars. Fortunately, the age fractions t/tms 
are much less sensitive to detailed calculations of stellar 
evolution. From the Ath and 5th block we see that the 
choice of tracks with rotation results in higher uncertain- 
ties on the fractions t/tms m the first evolutionary stages 
of the MS phase than the end of this phase. 



We can then conclude that mass determinations are 
not strongly sensitive to the type of evolutionary track 
used. On the other hand, the use of models for rotating 
stars that take into account the average rotational char- 
actersitics of fast rotators in the dwarf state of the MS 
leads to estimates of fractional ages which are not sensi- 
tive to the specific initial value of V Q around 300 km s _1 . 
This ensures that the models used in the present work lead 
to reliable statistical insights on global distributions and 
possible mass-dependencies of fractional ages at which the 
Be phenomenon occurs. It could be, however, suitable to 
proceed to more detailed iterations when discussing indi- 
vidual objects whose absolute ages are to be determined. 



4.2. Application to observed stars 

The entry parameters used to derive stellar masses and 
ages using relations ([[J and (J2J were considered with their 
la uncertainty bars: X — X ± a x {X Q stands for sur- 
face averaged logi or \ogg and logT e ff). Each interval 
(X —a x , X Q -\-a x ) was divided into 7 parts, so that the solu- 
tions of relations CJ-© and for each star interpolations 
in the HR diagrams were performed for all possible com- 
binations of individual sub-X^ entry parameters. Hence, 
for each star we obtained 8 2 solutions that determined 
the respectively r- and M-distributions of the solutions 
(most of them are not symmetrical). From these distri- 
butions we adopted the modes as the most probable re- 
sults, as well as the corresponding average la dispersion, 
to account for the related uncertainties. We note that the 
uncertainties affecting the apparent fundamental param- 
eters are those from the observed (Ai,D) quantities. We 
also have the uncertainty of (+0.06, —0.04) around the 
adopted rotation rate u = 0.88 that could affect the re- 
sults. Nevertheless, the global changes that will imply on 
the (t/tms, M/Mq) diagram the treatment of the funda- 
mental parameters with oj = or u> = 0.88, justify neglect- 
ing the small dispersion Slo = +0.06, —0.04. The apparent 
and surface average (T e s,logg) sets are given in Table 1. 
In this table we also reproduce the obtained ages, masses 
and MS age fractions derived using evolutionary models 
without and with rotation. 

Fig. Et) shows the HR diagram of the studied stars 
in terms of their apparent (logL,logT ff) parameters and 
where are sho wn also the evolutio nary tracks for non- 
rotating stars (Schaller et al 1 11 9921) . Fig. Eh) shows the 
HR diagram of the same stars, but in terms of their sur- 
face averaged (log L / L Q ,\ogT c g) quantities, where we as- 
sumed all stars rotate at w = 0.88. In this figure are also 
shown the evolutionary tracks for rotating objects that 
start evolving from the ZAMS as rigid rotators with equa- 
torial velocity V D = 300 km s^ 1 (MM2000). 

In spite of the widened MS phase of rotating stars 
and corrections made to their fundamental parameters 
for rotational effects, four program stars lie in the post- 
MS 'bright giant' region. They are HD23630, HD45910, 
HD183656 and HD217675, whose true Vsini in km s" 1 
and estimated inclinations i are (149;45°), (254;76°), 
(274;77°) and (272;90°) respectively. These objects are ref- 
erenced in the literature as binaries or multiple systems. 
Their apparently too 'evolved' character could be due to 
a merging effect of fundamental parameters from several 
components. 

The translation of the HR positions in Fig. E] into 
masses M/Mq and age fractions t/tms is given in Table 
1 and in Fig. EJ Fig. Et) shows the distribution of points 
(t/tmSj M/Mq) obtained for the program stars using the 
original or apparent fundamental pa rameters and the evo - 
lutionary tracks without rotation ijSchaller et alJll 992). 
The plotted error bars correspond to the la dispersions 
of the 8 2 t/tms- an d M/M©- solution distributions. Fig. 
Et>) shows the same type of distribution, but where param- 
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Fig. 5. HR diagrams of the studied Be stars, a) Apparent (log L/L®, logT c g) with evolutionary tracks without rotation, 
b) Surface averaged (logL/L0,logT e gf) and evolutionary tracks with initial an equatorial rotation velocity V = 300 
km s _1 for all masses. 



eters were obtained using surface averaged fundamental 
parameters for oj — 0.88 and models of stellar evolution 
with rotation calculated by MM2000 for V Q = 300 km s" 1 . 

We note that the uncertainty bars affecting the mass 
determinations are small, while those concerning the MS 
life fractions are in most cases quite large. This indicates 
that errors in the determined quantities due to approxi- 
mate choices of V can be concealed within the uncertain- 
ties related to the observations. 

In both diagrams of Fig. El points are spread over the 
whole interval of age fractions 0.0 t/tms Si 1-0, which 
suggest that the Be phenomenon may appear at any stage 
of the stellar evolution in the MS evolution phase. There 
is, however, a difference between solutions where rotation 
is taken into account and those where it is not. If we were 
not aware of rotational effects, Fig.Et) would suggest that 
86% of stars are above the t/tms = 0.5 limit. Fig. EJd) 
shows, however, that when the fast rotation of Be stars 
is taken into account, the fraction of stars in our sample 
above r ~ 0.5tms drops to 62%. 

Another important result appears when we separate 
the stars into massive (M > 12M Q ) and less massive 
ones (M ^ 12M ). We see then that in Fig.Eb) that the 
stars are distributed as follows: 

The mass-dependent division of Be stars regarding their 
position with respect to r = 0.5tms implies that the Be 



[upper 


MS half] 


30% 


65% 


[lower 


MS half] 


70% 


35% 






M ;> 12M 


M <> 12M 



'phenomenon in massive stars tends to appear on average 
at smaller t/tms a 9 e fractions than in the less massive 
stars. The same phenomenon is also suggested in Fig. [Ht>) 
by the stars with M 12M in the lower MS half region, 
as they all lie above a diagonal that starts among the more 
massive objects at t/tms — and ends at t/tms — 0.5 for 
M ~ 3M©. Since dwarf Be stars rot ate on average with 
V ~ 300 km s" 1 l|Zorec et al.ll2004|) . either we adopted 
V = 300 km s" 1 or V Q = f(M) ^ 300 km s" 1 . It can be 
shown that the initial angular velocity rate oj is higher 
the lower the stellar mass (cf. Sect. 13. 3|) . The fact that in 
the M <3 12M Q region stars are above this well defined 
mass-dependent slope might suggest that the mechanisms 
of angular momentum transport that accelerate the stellar 
surface up to a near/or critical velocity have relatively 
longer time scales as the stellar masse decreases. 

We also obtained masses and age ratios t/tms with 
the re-scaled evolutionary tracks for V a — f(M) 300 
km s _1 (cf. Sect. I3.4JI . The tilted character of the trend 
of points obtained is the same as the one presented in 
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Fig. 6. Age ratios t/tms of the studied Be stars against the mass, a) Parameters derived using evolutionary tracks 
without rotation; b) parameters derived using evolutionary tracks with rotation 



Fig-Eb)- There is, however, a slightly more pronounced 
downard drift of points. 

In Fig. we notice that 72% of the studied objects 
have masses from 3 to 12M . On the other hand, there is 
a striking lack of massive Be stars with M <; 16M Q . This 
lack can be due to several reasons: 1) smaller IMF for 
massive objects; 2) the fast evolution of massive stars, so 
that many of them in the solar neighbouring regions could 
have alreday migrated towards the supergiant phase; 3) 
the CE can be blown away by the radiatio n pressure of 
the 'potential' hottest Be stars <|Massalll97,4 



5. Discussion 

5.1. Comparison with previous works 

The HR diagram drawn in terms of the apparent BCD 
(Ai , D) parameters confir ms the tendency noti c ed by 
IJaschek et all (Il980f> and iHubert-Delplace et all l|l982|) 
that late-type Be stars can be on average slightly more lu- 
minous than early Be stars. This result can also be likened 
to the trend shown in Fig. Et) which concerns the appar- 
ent stellar parameters. 



Our finding of the Be phenomenon appearing at any 
ev olutionary phase is in a greement with s imilar assertions 
bv lMermilliodl l)l982|) and lSlettebakl l)l985|) based on stud- 
ies of Be stars in clusters. However, the scatter of points 
over the whole MS life span found by these authors can 
in part be due to the photometric data that were not cor- 
rected for CE perturbing effects and changes introduced 
by the fast rotatio n. 

IZorec fc Briol £l997) have concluded that the fre- 
quency of Be stars as a function of spectral type is roughly 
the same in all luminosity classes. These authors have 
dealt with a sample of field Be stars three times larger than 
the sample studied in the present work. However, they 
have gathered the studied objects in three wide groups 
of luminosity classes, where both evolutionary and rota- 
tional effects on the apparent luminosity class were mixed. 
Statistical corrections for CE flux excesses and for changes 
related to fast rotation were introduced. Since these cor- 
rections focused on the countings per group of luminosity 
class separately, they concerned only absolute magnitudes 
and spectral types. Moreover, in the number frequencies 
N(Be)/N(B+Be) against spectral type, not only do B stars 
without emission greatly outnumber Be stars, but they 
represent a highly heterogeneous group in physical charac- 
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teristics. Their evolution then has different characteristics 
to that of Be stars, which is thus a source of deviations in 
the count frequencies that c annot be ascrib e d to B e stars. 
Hence, the result found bv IZorec k, Briotl l)l997j) should 
be considered as intermediate to those shown in Fig. EJi) 
and Fig. Et>). 

When it comes to Be stars in clusters, statistics can 
be incomplete in the hotter and in the cooler extremes 
of spectral types. The scarcity of massive Be stars was 
already discussed in Sect. 14.21 To this add the possible 
exclusion of genuine massive Be stars in young clusters. 
In the accr etion pre-main-sequence (PMS) star formatio n 
paradigm l|Palla fc Stahlerlll993UBeech & Mitalaslll994h 
stars with masses M ^ 1 0M^ l(Behrend fc Maederll200lt 
iMaeder & Behrendl Eool may have a period of MS life 
when they can still be imbedded in the protostellar nebula. 
At these early evolutionary stages, there may be genuine 
fast rotators that have already formed, by mass-loss pro- 
cesses, some circumstellar disc. Thus, their observed emis- 
sion has two sources, the disc and the leftovers from the ac- 
cretion nebula. In spite of their mixed Be and Herbig AeBe 
character, these objects should not then be discarded from 
the Be statistics. For late type Be stars, their number is 
reduced because the low effective temperature maintains 
hydrogen atoms almost neutral, thus disabling possible 
existing CE to produce observable Balmer emission. Since 
the existing studies of frequencies of Be stars in clusters 
have not taken into account rotationally-induced effects 
in the spectral classifications, they should be compared 
with our results presented in Fig. EH). Thus, if in this fig- 
ure we disregard Be stars with masses M <; 12M and 
M <, 5Mq to mimic possible misdetections, we can see 
that the stars in the upper half of the diagrams widely 
outnumber those in the lowe r half. This readily accounts 
for the suggestion made by iFabregat fc Torre ionl |2000) 
that the Be phenomenon tends to appear in the second 
half of the MS life time. 



5.2. Tilted distribution of points 

Let us discuss briefly the tilted trend of points shown in 
Fig-Eb)- It is known that the mass-loss phenomenon in Be 
stars has two main components: a) winds with rates of the 
order of M ~ 10 _9 M Q yr _1 ; b) discrete mass-loss events 
of SM ~ 10~ 9 M(Ti underlying light outburs ts several times 
per year ((Hubert et al.ll20Q{l IZorecll2004h . Be stars with 
masses M <> 12M Q spend some 10 7 to 10 8 yr in the MS 
phase (MM2000), so they can be considered as evolving 
at a nearly constant mass regime. Although the mass-loss 
rates currently assumed for stars in the 12 < M/M & Ss 25 
mass interval are higher than quoted above, they cannot 
account for a total loss of the order of ~ 1OM during the 
MS phase to explain the slope of the upper edge of the 
strip of points in Fig. in terms of a drift towards the less 
massive side as the stars evolve from ZAMS to TAMS. 



The right-hand slope of the lower edge of the trend in 
Fig. Eh>) may suggest that B stars of low mass need to 



spend some mass-dependent time in the MS before they 
can display the Be phenomenon; i.e. the time needed for 
the surface velocity be spun up to f2/£! c ~ 0.9. 

We can also speculate on the tilted aspect of the dis- 
tribution of points in Fig. Eh) as produced by a mass- 
dependent distribution of initial equatorial velocities in 
the ZAMS, i.e. that massive stars start their MS phase 
with higher average rotational velocities than the less mas- 
sive ones relative to the respective critical velocities. While 
the suggestion of a mass- dependent evolution of surface 
rotational velocities needs a thorough theoretical study, 
some clues on the initial rotation in the ZAMS can be ob- 
tained from the basics of stellar structure and correlations 
drawn from observations. 

Evolutionary models of rotating stars 
l(Mevnet fc Maederl EoOfl l2002h suggest that we can di- 
vide the interior of stars in the ZAMS roughly into two 
regions, core and envelope, which rotate each at different, 
but respectively at near uniform angular velocities. The 
total stellar angular momentum J can then be written as: 



CO 

= k E M*R 2 n E x / 



k E (M* - M co ){R 2 - R 2 co )n E 



(8) 



where the k co and k E are the gyration radii; M co and M E 
are the masses; R co and R E the radii; tt co and tt E are the 
angular velocities of the core and envelope respectively; 
'*' stands for the whole star. To obtain from JHJ) an in- 
sight on the mass dependence of the initial velocity V , let 
us note that in MM2000 models we have M co /M* ^ 0.4, 
Rco/M* ^ 0.4, Q co /Q E ^ 3. Moreover, assuming that 
k co ~ 0.132 and k c ~ 0.04 as it would correspond to 'aver- 
age' polytropes of indices n = 1.5 for the convective core 
and n = 3.0 for the radiative envelope, we obtain that / 
in relation (|HJ) is fairly independent of mass and / ~ 1 . 

On the ot her hand, using observed Vs'mi values 
iKawalerl ljl987t ) has shown that on average the angular 
momentum of dwarfs depends on the stellar mass as: 



J ~ Mi 



(9) 



Since the initial MS radii of stars R scale with the stellar 
mass as: 



Ro ~ M 



0.6 



(M* Z 2M ), 



(10) 



from © for R = R Q the ratio of the average initial equato- 
rial velocity of stars V = fl E R normalized to the critical 
velocity V cr ~ (GM*/1.5i? ) 1 / 2 is related to the stellar 
mass as: 



fcrc 



(11) 



Knowing that for the whole range of stellar masses con- 
cerned by the Be phenomenon, models of the stellar inte- 
rior in the ZAMS are homologous, it follows that the gy- 
ration radius k E is independent of mass over the surveyed 
range of masses. Thus, <|1H gives us a first insight into 
the possible initial dependence of rotational velocities on 
mass, which implies that stars with masses A/* ^ 15M 
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can have initial rotational velocities 25% closer to the crit- 
ical one than those with masses M* Ss 5M Q . This mass 
dependence is small indeed. However, to account entirely 
for the noted tendency that hot stars reach the near/or 
critical rotation at smaller ratios t/tms than the cooler 
ones, there must still be some interaction between rota- 
tionally and evolutionary-induced changes on the stellar 
moment of inertia that can be felt as due to a small mass- 
related decrease of the gyration radius. 

Let us finally note that if we neglect the possible small 
mass-dependence of &e from M* ~ 15M Q to ~ 5M , re- 
lations IHJ and iJTlJ will imply that V r (15)/V r (5) ~ 1.6, 
while from observations of dwarf Be stars, which are on 
average at a later evolutionary stage than implied by V , 
we derive V(15)/V(5) ~ 1.3 (cf. Sect.lO)). This suggests 
that to obtain results of the kind shown in Fig. Eb), it 
would be better to use models with V = f(M) > 300 
km s" 1 (cf. Sect. [3~3Jl. rather than V Q = 300 km s^ 1 . 
Such a choice would imply the use of re-scaled evolution- 
ary tracks. We preferred, however, to keep the original 
models by MM2000, since they depict the effects of fast 
rotation in a more consistent way and because altogether 
they give reliable orders of magnitude of these effects. 

5.3. Comparison with model predictions 

The noticeable lack of Be stars with masses M 12M Q 
in the upper half of the MS can be explained as the nat- 
ural result of angular momentum loss produced by mass- 
loss, whose rate increases with the stellar mass. This loss 
reduces the surface equatorial rotational velocity conver- 
ting the star into a much lower rotating object (ui <§; 0.9). 
Conversely, the increase of the number of Be stars with 
masses M Ss 10M© in the upper half of the MS must be 
an effect of internal coupling, where the angular momen- 
tum is conveyed from the stellar core to the surface by the 
meridional circulation. The time scale of the meridian cir- 
culation is roughly T C j rc ~ tkh/v where tkh is the Kelvin- 
Hclmholtz time and rj is the ratio of the centrifugal force to 
the gravity. Since for Be stars it is r\ ~ 1, T C j rc ranges with 
the stellar mass as T circ ~ 4 x 10 6 (M/M Q ) -2 yr. It then 
becomes clear that the lapse of time to reach f2/f2 c ~ 1 is 
longer the smaller the stellar mass (MM2000 and Maeder 
2004, private communication), similar to model predic- 
tions. 

5.4. On the rotational velocities and star formation 
regions 

From hipparcos parallaxes we see that 84% of stars in our 
sample lie in a region within 500 pc of the Sun, of which 
62% are within 300 pc. It is also noteworthy that 45% of 
the more massive objects of this sample, hotter than 22000 
K, are within 300 pc and 35% are between 300 and 500 
pc. Only two stars are at d ~ 800 pc. The stellar sample 
studied should not then be characterized by strong dif- 
ferences in the initial metallicity. In fact, low metallicities 



migh t favor the fast rotation in some cases l|Maeder et alJ 
1999), but it cannot be the case for the massive Be stars 
of our sample. According to lMevnet fc Maederl (2002) low 
metallicity reduces the mass-loss rate, which favors the 
conservation of angular momentum in the stellar surface 
and s o, the existence of h i gher s urface rotational velocities. 
Also, iMaeder fe Mevnetl (2001) predict that the increase 
of the O / O c ratio towards 1 in massive stars is faster in 
models with low metallicity. It would then be interesting 
to test this prediction by obtaining diagrams like those in 
Fig. for Be stars in environments with quite different 
metal abundance, an aim that we will pursue in subse- 
quent work. 



5.5. On the initial conditions 

The calculated properties of an evolving object with rota- 
tion in the MS depend on the assumed initial conditions. 
The results given in the present work depend on predic- 
tions made for stars which began evolving in the ZAMS as 
rigid rotators. This choice may, however, be not the only 
possible. Similarly to ZAMS, which is likely a computa- 
tional 'landmark' l|Endal fc Sofialll98lh . rigid rotation in 
the ZAMS may be a simplification too. On the one hand, 
the necessity and/or definition of a ZAMS for massive 
stars is not clear, as it happens in the accretion parad igm 
of star formation jPalla fc Stahlerlll993tlBeech fc Mitalad 
Il994t IMaeder fc Behrendll2002|) . On the other hand, due 
to hydrodynamical instabilities, an initial rigid rotation 
switc hes rapidly (some 10 4 y r ) into a differential rota - 
tion l|Denissenkov et alJll99ft iMevnet fc Maederl l2000h . 
Such a differential rotation may then be present before 
the 'ZAMS' phase. 

In the classical PMS evolution frame, based on the 
contraction of a constant mass sphere until the gravita- 
tional energy release increases the central temperature 
enough to trigger the nuclear reactions, rigid rotation 
in the ZAMS was generally justified because: a) dynam- 
ical stability against axisymmetric per turbations coul d 
be warranted for rigid stellar rotators !|Fuiimotd ri987tk 
b) it was assumed that in the PMS the full convection 
phase stars become rigid rotators. However, recent 2-D 
hydrodynamical calculations show that convection does 
not maintain rigid rotation, but it rather produces an in- 
ternal angular velocity distribution profile £l(m) cx vo~ v 
(zu = distance to the rotation axis), which is interme- 
diate between complete redistribution of specifi c angular 
momentum (p = 2) and rigid rotation (p = 0) l|Deupred 
I1998L l200fi 1200 Furthermore, in the accretion forma- 
tion scheme th e star is a mass and angular momentum 
gaining object Jpalla fc Stahlerlll993T : lMercer-Smith et alJ 
1999; iTerebev et all l984). Meanwhile a non-rotating star 
gains mass, it first undergoes a full-convection period, then 
a radiative core is developed, the star swells until a tem- 
porary full radiative state is attained and finally the core 
becomes convective. During this period the forming star 
stretches and contracts, but it can also undergo magnetic 
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interaction with the accretion disc llSteDienll200i . These 
phenomena can produce an uneven distribution of the an- 
gular momentum inside the star. There is also some de- 
coupling of stellar internal regions, since the time scales of 
angular momentum redis tribution in convectiv e and radia- 
tive zones are different (lEndal fc Sofialll98ll) . If the star 
acquires mass through an accretion disc, which is prob- 
ably in Ke plerian rotati on, a huge gain of angular must 
take place (|PacketJll98l : 

j(t) w [GR{t)M{i)\ x / 2 M(t) (12) 

where G is the gravitional constant, R and M are the 
time-dependent radius and mass of the star respectively 
and M is the mass accr etion rate \M ~ 10~ 5 (M /M^)^ 
M Q /yr with <f> ~ 1.0-1.5 l|Maeder fe Behrendl2002|) ] . Since 
the rigid rotation implies a very low amount of rotational 
energy, it cannot be then excluded that a considerable 
load of angular momentum may exist somewhere deep in- 
side the star, so that an a priori assumption of rigid rota- 
tion can be d ifficult to justify. However, it was noticed by 
ISpruitl l|l999ll2002]) that a powerful enough differential ro- 
tation can create magnetic fields. Different regions inside 
the star could then be "locked" to reco ver some degree of 
rig id rotati o n llMaeder fe Mevnedl2003|) . 

IStepiehl l)2002|) has suggested that given an appropri- 
ate range of surface magnetic fields, stars may gain angular 
momentum through an effe ctive magnetic accretion in the 
PMS nha,se. lSter)ie! notes that the rotation can be 

faster the more massive the star, because the interaction 
with the circumstellar matter lasts less time. Nevertheless, 
the final balance between losses and gains of angular mo- 
mentum produced by interactions with the circumstellar 
environment were not definitively established, nor were 
its consequences on the internal rotation law of the star. 
In particular, very little was said about the amount of 
rotational kinetic energy the star is left with after these 
interactions. 

Rigid rotation puts an upper limit onto the amount 
of rotational kinetic energy /C a star can store. At rigid 
critical rotation, an early-type star has £ c — /C C /|W| ~ 
0.015 (W = gravitational potential energy). In a star with 
internal differential rotation, the same surface rotations 
may correspond to higher values of £, which may then 
carry stronger stellar deformations, gravity darkening ef- 
fects and internal hydrodynamical instabilities. For an or- 
der of magnitude estimate, Table0]gives the rotational ki- 
netic energy ratios /C(p)//C(p=0) and energy ratios £ (p) = 
IC(p)/\W(p)\ for an internal rotational law H(vj) oc m~ p , 
assuming that the stellar surface rotates at f2 s /fi c = 0.9 
as occurs on average for Be stars. These values were ob- 
tained using two d imentional models of stellar structure 
I Zorec et al.lll988a| . We see that for a mild differential 
rotation p = 0.4 it is £ > £ c and that for p — 0.7, an 
average rotation law set by convection, £ is nearly twice 
as high as for a critical rigid rotator. In such a case there 
is a lowering of the core bolometric luminosity that ranges 
from roughly 17% at masses M ~ 3OM to 27% in masses 



Table 4. Kinetic energy ratios and £ ratios calculated for 
different values of p 



p 


/C(p)//C(p=0) 


£(p) 


0.4 


1.52 


0.018 


0.7 


2.17 


0.026 


1.0 


3.26 


0.038 



M ~ 3M dClemendll97l . A much complicate relation 
must then exist between spectra, masses and stellar ages 
than treated in the present work llCollins fc Smith! f l 985: 
IZoredll986llZorec et al.lll987Lll988blaLll990tlZoredll992() . 

6. Conclusions 

In this paper we have studied a sample of 97 field Be 
stars, most of which are at distances d < 500 pc from the 
Sun, so that they can be considered more or less homo- 
geneous regarding their initial metallicity. All these stars 
were observed in the BCD spectrophotometry system to 
have photospheric spectral signatures as much as possi- 
ble free of CE emission/absorption perturbations. The ap- 
parent fundamental parameters derived from the observed 
BCD (Ai,D) quantities, i.e. parameters reflecting the av- 
erage rotationally-perturbed photosphere shown by the 
projected stellar hemisphere towards the observer, were 
translated into pnrc and averaged fundamental parame- 
ters, "pnrc" is the acronym for parent non-rotating coun- 
terparts, or parameters that correspond to homologous 
non-rotating stars. The averaged fundamental parame- 
ters correspond to averages over the whole stellar surface. 
We have assumed that the studied B e stars rotate with a n 
angular velocity ratio 0/f2 c = 0.88 ijFremat et alJl2005^ . 
Note the difference between averaged and apparent, the 
last representing a sort of average spectrum emitted by 
the 'observed' stellar hemisphere. The averaged parame- 
ters are the only quantities that can be used to interpolate 
stellar masses and ages in the evolutionary tracks. 

The present contribution represents one of the first at- 
tempts to derive stellar masses and ages of Be stars by 
using simultaneously model atmospheres and evolution- 
ary tracks both calculated for rotating objects. According 
to the statistical average of true rotational velocities V 
of dwarf Be stars, the evolutionary models used are for 
ZAMS equatorial rotational velocity V — 300 km s _1 in 
all masses. For all stars we derived the mass and stellar 
ages r normalized to the respective time that each rotat- 
ing star can spend in the main sequence phase tms- As 
a consequence of effects of the rapid rotation described 
by the models used, we obtained a trend of points in the 
(t / tms , M / Mq) diagram, which implies that: 

a) there are Be stars spread over the whole age interval 
^ t/tms 1 in the main sequence evolutionary phase; 
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b) in massive stars the Be phenomenon tends to be 
present at lower t/tms age ratios than in the less massive 
stars. 

The lack of massive stars in the upper MS can be due 
to loss of angular momentum through mass-loss which pro- 
duces a strong decrease of the stellar surface f2/fi c ratio. 
On the other hand, the increase of the number of less mas- 
sive Be stars in the upper MS can be explained in terms of 
angular momentum transport from the core to the surface 
carried by the meridional circulation, which has longer 
time scales the lower the stellar mass. 

Arguments based on the distribution of the total an- 
gular momentum of dwarf stars against mass reveal that 
the massive stars may start evolving from the ZAMS with 
a slightly higher V o /V than the less massive ones. 
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Table 5. Program stars, observed and derived parameters 
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Table 6. Program stars, observed and derived parameters 
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